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NUCLEAR STAR-FORMING RING OF THE MILKY WAY: SIMULATIONS 
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ABSTRACT 

We present hydrodynamic simulations of gas clouds in the central kpc region of the Milky Way 
that is modeled with a three-dimensional bar potential. Our simulations consider realistic gas cooling 
and heating, star formation, and supernova feedback. A ring of dense gas clouds forms as a result of 
X1-X2 orbit transfer, and our potential model results in a ring radius of ~ 200 pc, which coincides 
with the extraordinary reservoir of dense molecular clouds in the inner bulge, the Central Molecular 
Zone (CMZ). The gas clouds accumulated in the CMZ can reach high enough densities to form stars, 
and with an appropriate choice of simulation parameters, we successfully reproduce the observed gas 
mass and the star formation rate (SFR) in the CMZ, ~ 2x 10 7 Mq and ~ 0.1 M® yr . Star formation 
in our simulations takes place mostly in the outermost X 2 orbits, and the SFR per unit surface area 
outside the CMZ is much lower. These facts suggest that the inner Galactic bulge may harbor a mild 
version of the nuclear star-forming rings seen in some external disk galaxies. Furthermore, from the 
relatively small size of the Milky Way's nuclear bulge, which is thought to be a result of sustained 
star formation in the CMZ, we infer that the Galactic inner bulge probably had a shallower density 
profile or stronger bar elongation in the past. 

Subject headings: Galaxy: center — Galaxy: nucleus — galaxies: star formation — galaxies: ISM — 
ISM: kinematics and dynamics 



1. INTRODUCTION 

CO emission surveys along the Galactic plane reveal 
that molecular gas is abundant in the plane down to 
a Galactocentric radius r of ~ 3 kpc (see e.g., Dame, 
Hartmann, & Thaddeus 2001). Inside this radius, the 
amount of molecular emission greatly reduces, as do star 
formation activities. However, a significant amount of 
molecular g well as various evidences of recent star 
formation, appears again inside r of ~ 200 pc (Mor- 
ris & Serabyn 1996, hereafter MS96; Yusef-Zadeh et al. 
2009, hereafter YZ09). This region of abundant molec- 
ular gas is called the Central Molecular Zone (CMZ). 
The total gas mass inside r — 250 pc is estimated to be 
1.7-5 x 10 7 M Q (Dahmen et al. 1998; Pierce-Price et al. 
2000; Launhardt, Zylka, & Mezger 2002). 

The gas in the CMZ must have migrated inward 
from the Galactic disk through interactions with other 
clouds, field stars, and magnetic fields (MS96). A mixed 
atomic/molecular layer inside r ~ 3 kpc turns into 
a mostly molecular, high-density medium at the CMZ 
(Liszt & Burton 1978; MS96). Binney et al. (1991) 
hypothesize that this abrupt change in gas content and 
density is caused by a transition between so-called X\ 
and X2 orbits near r = 200 pc. X\ is a family of sta- 
ble, closed orbits elongated along the bar's major axis, 
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whereas X2 orbits are elongated along the bar's minor 
axis, much deeper in the potential (Contopoulos & Pa- 
payannopoulos 1980). The spiraling-in gas clouds would 
generally follow the X\ orbits until the orbits become 
sharply cusped and even self-intersecting at the inner- 
most X\ region. Gas clouds experience compression and 
shocks near the cusps of these orbits, lose angular mo- 
mentum, and plunge to the X2 orbits inside. The in- 
falling gas clouds can additionally lose angular momen- 
tum when they collide with those already on the X2 
orbits. This compression and subsequent cooling will 
transform the gas into molecular form, and the molecular 
clouds accumulated in the X2 region would correspond 
to the CMZ. 

The abundance and high densities of molecular gas in 
the CMZ naturally lead to star formation. Indeed the 
CMZ harbors two extraordinary, young (few Myr), mas- 
sive (a few 10 4 M Q ) star clusters, the Arches and the 
Quintuplet (Figer et al. 1999; Kim et al. 2006, among 
others) as well as a population of young stellar objects 
(YZ09). Serabyn & Morris (1996) propose that star for- 
mation has been occurring in the CMZ throughout the 
lifetime of the Galaxy, and the resulting stellar popula- 
tion is evident as the central r~ 2 cluster or the nuclear 
bulge, a prominent cusp in the stellar density profile, 
whose extent is similar to that of the CMZ. 

There have been a few numerical studies on the hydro- 
dynamics of the gas layer in the Galactic inner bulge. 
Jenkins & Binney (1994) examined the hypothesis of 
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Binney et al. (1991) with two-dimensional (2-D) sticky 
particle simulations and confirmed that the transition 
of gas motion from X\ to X 2 orbits indeed takes place 
in a bar potential. Lee et al. (1999) and Englmaier 
& Gerhard (1999) were able to observe the same phe- 
nomenon in their 2-D smoothed particle hydrodynam- 
ics (SPH) simulations. Rodriguez-Fernandez & Combes 
(2008) performed 2-D hydrodynamic simulations (parti- 
cle motions followed in 3-D) with a more sophisticated 
potential model for the inner Galaxy in an attempt to 
reproduce the observed velocity structure of the CMZ. 
All these studies primarily concentrated on the behavior 
of the gas flows, and were not able to follow the density 
and temperature evolution of the gas clouds and subse- 
quent star formation^ due to the limits in their numeri- 
cal schemes and/or computing power at those times. 

In the present Letter, we present 3-D hydrodynamic 
simulations for the motions of gas clouds and star forma- 
tion in the inner bulge of the Galaxy. Our simulations 
consider self-gravity, gas cooling and heating, and super- 
nova (SN) feedback. We show that a bar potential can 
indeed compress and cool the gas sufficiently enough to 
form stars and that the obtained star formation rates 
(SFRs) are consistent with observations. 

2. METHODS AND MODELS 

We use a parallel A-body/SPH co de ASURA (Saitoh 
et al. 2008, 2009) for our simulations. We use an open- 
ing angle of 9 = 0.5 for a cell opening and a softening 
length of 3 pc, and the kernel size of an SPH particle 
is determined by imposing the number of neighbors to 
be 32 ± 2. We implement a cooling function by Spaans 
& Norman (1997) for gas with a solar metallicity for 
a temperature range of 10-10 8 K, and uniform heat- 
ing from far-ultraviolet (FUV) radiation is considered. 
The standard radiation strength used in our simulations 
is 100 Go0 which is ~ 60 times the solar neighborhood 
value (the SFR per unit surface area in the CMZ is ~ 60 
times larger than that in the solar neighborhood; YZ09). 
The effect of SN feedback is implemented by supplying 
thermal energy into the surrounding 32 gas particles. 

A star particle is spawned when a gas particle satisfies 
all of the following four conditions: 1) hydrogen number 
density is larger than the threshold density, n t h, 2) tem- 
perature is lower than the threshold temperature, Tin, 
3) the flow is converging (V • v < 0), and 4) there is 
no heat supply from nearby SN explosions. We adopt 
n t h = 100 cm~ 3 and T t h = 100 K as our standard 
threshold values following Saitoh et al. (2008, 2009), 
who showed that only such a high density and low tem- 
perature threshold can reproduce the structural charac- 
teristics of galactic disks. The local SFR is assumed to 
be proportional to the local gas density divided by the 
local dynamical time, 

OX ^dyn 

where C* is the star formation efficiency (SFE) param- 
eter, and we adopt C* = 0.033 as our standard value 

2 A recent Milky Way simulation by Baba, Saitoh, &; Wada 
(2010), whose main interest was gas kinematics, does follow ther- 
mal evolution and star formation of the gas. 

3 G = 1.6xl0" 3 ergcm- 2 s- 1 is the Habing field (Habing 1968). 



Table 1 

Simulation Parameters and Results 



Sim 


M a,i 
(1O 7 M ) 


C, 


^FUV 

(Go) 


M a,X2 
(ld f M ) 


SFR X2 
(Moyr- 1 ) 


1 


10 


0.033 


100 


1.33 


0.060 


2 


20 


0.033 


100 


2.03 


0.116 


3 


10 


0.0033 


100 


1.85 


0.045 


4 


10 


0.033 


1000 


1.11 


0.057 



Note. — M g i is the initial total gas mass, Jftjv is the FUV ra- 
diation strength, and M„ t x2 and SFRx2 are the gas mass and the 
star formation rate in the X2 region averaged for T = 200—600 Myr. 
The initial number of gas particles is 200,000 for simulations 1, 3, 
& 4, and 400,000 for simulation 2. 



following Saitoh et al. A newly spawned star particle is 
set to have one-third of the original gas particle mass. 
Detailed discussion on the spawning of star particles and 
SN feedback is given in Saitoh et al. 

Following Zhao, Spergel, & Rich (1994), we adopt an 
m = 2 bar with a power-law density profile for the inner 
Galactic bulge: 

p = p°(^) [i-nM)L (2) 

where 

Y(6, <f>) = -620^20 (cos 9) + b 22 P 22 (cos 0) cos 20, (3) 

and P is the associated Legendre function. These are 
linear combinations of spherical harmonic functions of 
the I = 2, m = 0,2 modes. &20 determines the degree of 
oblateness/prolateness while &22 determines the degree of 
non-axisymmetry. Our standard simulation (simulation 
1) has the following parameters: a — 1.75, b 20 = 0.3, 
622 = 0.1, po = 40 Mq pc~ 3 , r = 100 pc, and a bar 
pattern speed of 63km/s/kpc. This set of parameters 
gives a bar axis ratio of 1:0.74:0.65 for the isodcnsity 
surface that intersects points [x = 0, y = ±200 pc, z = 
0], an X\-X% transition at r ~ 200 pc, inner Lindblad 
resonance at 660 pc, and a corotation radius of 2.8 kpc. 
Enclosed masses inside 200 pc and 1000 pc are 10 9 M© 
and 7 x 10 9 M Q , respectively. 

Our simulations initially have gas particles only, which 
are distributed on X\ orbits whose semi-minor axes range 
from 300 to 1200 pc. The number of particles on each 
orbit is proportional to the length of the orbit, and on a 
given orbit, particles are distributed over the same time- 
interval. The initial vertical distribution of the gas par- 
ticles follows a Gaussian function with a scale height of 
40 pc. Our standard simulation has 2 x 10 5 gas particles 
and a total gas mass of 10 8 M Q , thus each gas particle 
initially has a mass of 500 M Q (the latter is true for all 
of simulations, thus star particles in our simulations typ- 
ically have ~ 170 Mq). The initial velocities in the plane 
follow the motions on the X\ orbits, while the vertical 
component is initially set to zero. We evolve such initial 
setup without cooling, heating, star formation, and SN 
feedback for the first 50 Myr to obtain a hydrodynami- 
cally relaxed particle distribution. 

3. RESULTS 

Three snapshots of our standard simulation are pre- 
sented in Figure [Q for T = 0, 150 and 500 Myr. The 
sense of galactic rotation in the snapshots is clockwise. 
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Figure 1. Surface density maps of gas and (accumulated) star particles in simulation 1 at T = 0, 150, and 500 Myr. Length units are in 
parsecs, and the gray scale represents densities linearly from 3 to 300 Mq pc -2 . The major axis of the bar potential is aligned along the 
y-axis. 
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2. Evolution of total (gas and stars), gas, and stellar masses (panels a through c), and the star formation rate (panel d) in the Xi 
lines) and X2 (solid lines) regions of simulation 1. After the initial relaxation period (T > 200 Myr), most of star formation takes 



Figure 

(dashed 

place in the X2 region (CMZ), and the obtained SFR in the CMZ is consistent with the observed values by YZ09, 0.01-0.1 Mq yr 



The middle panel (T = 150 Myr) shows that gas parti- 
cles near [x = 0, y = ±1250 pc] quickly undergo shocks 
and compressed clouds plunge to the X2 orbits along the 
"dust lanes" in the leading inner edges of the X\ orbits. 
The gas accumulated on the X 2 orbits, i.e., the CMZ, 
initially forms an elongated ring of clouds at r ~ 200 pc. 
The masses of gas clumps and streams in the CMZ of our 
simulations, 10 5 -10 6 M Q , agree with those estimated for 
some clumps from radio observations (e.g., Stark et al. 
1991). Gas density is the highest in this ring, and this is 
where star formation is the most active. 

The right panel (T = 500 Myr) shows that as the gas 
infall and star formation in the nuclear ring continue, 
the radial distribution of gas and stars in the X2 re- 
gion broadens (mostly inwarcQ). However, star forma- 
tion still takes place mostly near r ~ 200 pc, because 
newly infalling gas from the dust lanes collides the ex- 
isting nuclear ring of gas at its outer rim. We note that 
only ~ 15 % of the gas flowing down the dust lanes di- 
rectly enters the CMZ, while the rest flows back to the 
X\ region. This fraction is consistent with that found in 
Regan, Vogel, & Teuben (1997; 10-25 %). 

Evolutions of gas, stellar, and total masses on X\ and 

4 The gas particles in the CMZ lose orbital energy via shear 
viscosity from differential rotation and gradually migrate inward. 
Some particles should gain angular momentum and the CMZ 
should broadens outward as well, but those particles would merge 
with the part of the dust lane stream that returns back to the X\ 
region. 



X2 orbits are shown in Figure O The mass transfer from 
X\ to X2 regions takes place quite rapidly between T ~ 
100 and ~ 200 Myr and then goes into a slower, stable 
phase (panel a). While the total mass in the X2 region 
grows nearly constantly after T = 200 Myr, the gas mass 
in the X2 region, M 9t x2, is held at a nearly constant 
value and then decreases very slowly (panel b). This 
implies that star formation in the X2 region becomes 
so efficient that all of the newly supplied gas turns into 
stars once the gas mass in the X2 region reaches a certain 
critical value (~ 1.5 x 1O 7 M ) at ~ 200 Myr. 

Both the Xi and X 2 regions initially (T < 200 Myr) 
have high SFRs (see panel d), but this is because when 
the star formation takes place for the first time (T = 
50 Myr for the X\ region, ~ 100 Myr for the X2 region), 
there is no or not enough SN feedback that regulates star 
formation. Thus this initial relaxation period should be 
disregarded if one assumes a continuous gas supply from 
the outer region (an intermittent gas supply may result 
in a repetition of a brief, high-SFR phase followed by 
a longer, lower-SFR phase). While the SFR in the X\ 
region continuously decreases as the amount of gas there 
becomes depleted, the SFR in the X2 region becomes 
quite stationary after the initial relaxation period, and 
ranges between 0.04 and O.O9M yr _1 . This range is 
in excellent agreement with recent SFRs estimated for 
the CMZ from mid-infrared observations by YZ09, 0.01- 
0.1 M s yr" 1 . 

Simulation 1 clearly shows that 1) under a bar po- 
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tential modeled for the inner bulge of the Milky Way, 
gas inflow along the Galactic plane leads to the forma- 
tion of a gas ring, whose radial extent is consistent with 
the CMZ, and 2) the gas ring has high enough densities 
to form stars at a rate consistent with observed values. 
The parameters of our potential model (po, ro, a, &20 , 
and 622) were chosen such that the transition between 
X\ and X 2 orbits appears at ~ 200 pc, therefore the ex- 
tent of the gas ring in our simulation is a result of the 
potential parameters that we chose. 

On the other hand, it is reasonable to expect that the 
SFR in the CMZ is determined by the gas influx rate 
from the X\ to X2 regions, and that this influx rate is 
dependent on the initial gas surface density (or initial 
gas mass) in the X\ region. The fact that the SFR from 
simulation 1 is consistent with the observed values im- 
plies that the adopted initial gas mass in our X\ region 
(major axes from 300 to 1200 pc) properly reproduces 
the actual gas influx rate from the X\ to X2 regions in 
the GC. 

However, the average M g> x2 after relaxation (200 < 
T < 600 Myr) in simulation 1, 1.3 x 10 7 M Q , is 20 % to 
four times smaller than the estimated M g _x2 values from 
various observations, 1.7-5 x 10 7 M© (the largest mass 
estimate is derived from dust emission and is sensitive 
to the assumed dust temperature in the GC, while the 
smallest estimate is obtained from C 18 molecular line 
emission; see Section [T] for references). 

To see the dependence of M 9t x2 on initial M g< xi, we 
performed simulation 2, which starts with a twice as large 
initial M 9r xi than simulation 1. Simulation 2 indeed 
results in ~ 2 times larger SFRs than simulation 1 as 
expected, and ~ 1.5 times larger M 9: x2 (see Table [l}. 
Now, not only the SFR but also M 9t x2 is consistent with 
observations. Note that the CMZ mass is not as sensitive 
to the gas influx rate as the SFR is. 

Simulation 3, whose SFE parameter (G*) is 10 times 
smaller than simulation 1 in consideration of the ex- 
treme star formation environment in the GC such as 
strong tidal forces/magnetic fields and elevated gas tem- 
peratures, shows another way to increase M 9} x2 — it 
gives ~ 25 % smaller SFR, but ~ 40 % larger M ffiX 2 
(1.9 x 10 7 Mq) than simulation 1. The increase in M 9t x2 
is due to the fact that the star formation in the CMZ 
now requires higher gas densities than in simulation 1 to 
compensate for the lowered G* in Equation ((TJ. From 
our experiments with simulations 1 through 3, we con- 
clude that one can fit the observed gas mass and SFR 
in the CMZ simultaneously by controlling the gas influx 
rate from the outer bulge (i.e., initial M 9i xi in our sim- 
ulations) and the SFE (G*). 

The actual SFR in the whole X\ region of the Milky 
Way is not known yet, as no systematic observations have 
been made for the region. Simulations 1 through 3 show 
some star formation in the X\ region, but it is not possi- 
ble to accurately calculate steady-state SFRs in the X\ 
region from these simulations, because our simulations 
do not have continuous gas supply from the outer bulge. 
The total SFRs in the X\ region of simulations 1 through 
3 are only fractions of those in the X 2 region, implying 
that the SFR per unit volume is much higher in the X 2 
region (the volume occupied by the X\ orbits is ~ 30 
times larger than that occupied by the X 2 orbits in these 
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Figure 3. Evolution of gas and stellar masses in the Xi (dashed 
lines) and X2 (solid lines) regions of simulation 4. The increased 
Go effectively suppresses star formation in the X\ region, where 
the gas density is lower. 

simulations) . 

In many nuclear starburst galaxies, star formation 
takes the form of a ring, i.e., a nuclear star- forming ring. 
If the Milky Way is such a case, star formation in the in- 
ner bulge would be confined in the CMZ and there should 
be no or very little star formation taking place outside 
the CMZ. We find from simulation 4 that one way to 
achieve this is to impose a higher FUV radiation field. 
Simulation 4, whose Go is 10 times larger than simula- 
tion 1, gives nearly the same gas mass and SFR in the 
X 2 region as simulation 1, but almost entirely suppresses 
star formation in the X\ region^ (see Figure [3]) . It is not 
certain whether the radiation strength in the inner bulge 
can actually reach as high as 1000 Go, but other physical 
mechanisms that can quench star formation in the GC, 
but are not considered in our simulations, such as strong 
tidal forces and magnetic fields, may also contribute to 
the suppression of star formation in the X\ region as seen 
in simulation 4. 

Our star formation and SN feedback models success- 
fully reproduce the Kennicutt-Schmidt (K-S) relation. 
At T = 400 Myr, the CMZ in simulation 1 has a gas 
surface density of 100M©/pc 2 and a surface SFR of 
0.45 Mq yr _1 kpc 2 , the latter of which is 0.45 dex higher 
than the mean K-S relation, whereas the CMZ in simu- 
lation 3 has a surface SFR that is 0.12 dex higher than 
the mean K-S relation. Considering a rather large scat- 
ter in observed SFRs around a mean K-S relation (rms 
scatter ~ 0.3 dex; Kennicutt 1998), the SFRs from our 
simulations are consistent with the K-S relation. 

4. THE CMZ AND THE NUCLEAR BULGE 

If the star formation in the inner bulge of the Milky 
Way is indeed confined in the CMZ, and if it has been 
sustained for the lifetime of the Galaxy as suggested 
by Figer et al. (2004), the resulting stellar population 
should manifest in a certain way. Indeed, the infrared 
intensity profiles in the GC have revealed a prominent, 
central stellar cusp with an outer extent of about 0.7 deg 
(~ 100 pc) on top of the shallower profile made by the 
Galactic bulge (Serabyn & Morris 1996, and references 
therein). Serabyn & Morris propose that this 'nuclear 
bulge' is the resulting stellar population from the sus- 

5 The suppression of star formation by FUV radiation is much 
more effective in less dense gas. 
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tained star formation in the CMZ0 

The extent of the nuclear bulge is somewhat smaller 
than that of the CMZ, and this may indicate that the 
radius of the CMZ was smaller in the past. From our 
further simulations with various density slopes (a) and 
bar elogations (622), we find that a smaller CMZ radius 
can be obtained when the density profile is shallower 
(a < 1.6; the same amount of orbital energy loss at the 
cusps of X\ orbits will plunge the gas further inside), or 
when the bar is more elongated (&22 0.14; the gas at 
the cusps of X\ orbits will lose larger orbital energy). 
Therefore, the relatively smaller size of the nuclear bulge 
may be an important clue on the secular evolution of the 
Galactic center. 

Some of the accumulated gas in the CMZ will keep mi- 
grating inward and reach the CircumNuclear Disk (CND) 
of molecular clouds at r — 2-5 pc, and a part of it will 
sink further to form stars in the central parsec or feed 
the central supermassive black hole. One possible way 
to transfer gas from the CMZ to the CND would be a 
'nested bar,' as examined by Namekata et al. (2009). 
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